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Abstract 

In this paper we focus on the water waves problem for uneven bottoms on a two-dimensionnal 
domain. Starting from the symmetric Boussinesq systems derived in [Chazel, Influence of bot- 
tom topography on long water waves, 2007], we recover the uncoupled Korteweg-de Vries (KdV) 
approximation justified in [Schneider- Wayne, The long-wave limit for the water-wave problem. 
I. The case of zero surface tension, 2002] for flat bottoms, and in [Iguchi, A long wave ap- 
proximation for capillary- gravity waves and an effect of the bottom, 2007] in the context of 
bottoms tending to zero at infinity at a substantial rate. The goal of this paper is to investi- 
gate the validity of this approximation for more general bathymetries. We exhibit two kinds of 
topography for which this approximation diverges from the Boussinesq solutions. A topograph- 
ically modified KdV approximation is then proposed to deal with such bathymetries, where 
topography-dependent terms are added to the solutions of the KdV equations. Finally, all the 
models involved are numerically computed and compared. 

Key words: Water waves, free surface flows, uneven bottoms, bottom topography, long waves, 
Korteweg-de Vries approximation, Boussinesq models. 



Introduction 

The water waves problem consists in describing the evolution of the free surface and 
velocity field of a layer of ideal, incompressible and irrotationnal fluid under the only in- 
fluence of gravity. The governing equations - also called free surface Euler equations - are 
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fully non-linear and non-strictly hyperbolic, and their direct study and computation re- 
mains a real obstacle. Many authors such as Nalimov ([44], 1974), Yoshihara ([57], 1982), 
Craig (pi], 1985), Wu ([55], 1997 and [56], 1999), Ambrose-Masmoudi (0], 2005) and 
Lannes ([38], 2005 and [2], 2007) have successfully tackled the problem of well-posedness 
of these equations. Nevertheless, the numerical computation of these solutions remains 
a tough task, especially in 3-D - see the works of Grilli et al. ([25], 2001) and Fochesato- 
Dias ( 23 , 2001). 

An alternative way to describe these solutions and their time behaviour is to look for ap- 
proximations via the use of asymptotic models. Such models are usually derived formally 
from the water waves problem by introducing dimensionless parameters. Making some 
hypothesis on these parameters reduces the framework to more limited physical regimes 
but allows the construction of asymptotic models. In this work, we focus on the so-called 
long waves regime. In this regime, the ratios e = a /ho and /x = h^/X 2 where a denotes 
the typical amplitude of the waves, ho the mean depth and A the typical wavelength, are 
small and of the same order. Many models can be found in the litterature corresponding 
to this regime. Among them, we can quote the works of Boussinesq ([E], 1871 and [15], 
1872) who was the first to propose a model that take into account both nonlinear and dis- 
persive effects, the unidirectional models such as the Korteweg-de Vries (KdV) equation 
([36], 1895), the Kadomtsev-Petviashvili (KP) equation ([31], 1970) and the Benjamin- 
Bona-Mahony one ([5], 1972). These historical models have been considerably studied 
and generalized, and their justification has been investigated among others by Craig 
(PI], 1985), Schneider- Wayne ([50], 2000), Bona-Colin-Lannes ([12], 2005), Lannes-Saut 
(gD], 2006) for flat bottoms, and Iguchi ([2B], 2006), Chazel ([M], 2007), Alvarez-Lannes 
([2], 2007) for uneven bottoms. 

If we focus more specifically on the KdV approximation, we emphasize that the investiga- 
tion of this model on uneven bottoms is not recent : this subject has been tackled among 
others by Ostrovskii-Pelinovskii ([48], 1970), Kakutani ([32], 1971), Johnson ( 28J, 1973), 
Miles (.43], 1979) and Newell ([H], 1985). However, their approach is quite different com- 
pared to the one chosen here, as will be clarified later in this paper. The more recent 
articles of Iguchi ( 26 , 2006) and Chazel ([16], 2007) are somehow the starting point of 
the present work. In [26 , Iguchi derived two different models : a coupled KdV system for 
relatively general bottom topographies, and a uncoupled KdV system for bathymctries 
decaying at a substantial rate at infinity. In [16] , the author derived two classes of sym- 
metric Boussinesq systems for two bottom topography scales, corresponding respectively 
to slightly and largely varying bottoms. The aim of this work is to recover the uncoupled 
KdV model - also justified by Schneider- Wayne [50] for flat bottoms - starting from any 
previous Boussinesq system proposed in [16j for slightly varying bottoms, to discuss its 
validity regarding the bottom topography, and to propose an alternative. 



Formulation of the problem 

In this paper, we work in two dimensions : x corresponds to the horizontal coordinate 
and y to the vertical one. We denote by (t, x) — > rj(t, x) and x — > b(x) the parametrizations 
of the free surface and bottom, defined respectively over the surface y = and the mean 
depth y = —ho at the steady state. The time-dependant domain ft t of the fluid is thus 
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Figure 0.1. Representation of the fluid domain 



y =Ti(t,x) 



n, 




= -h +b(x) 



taken of the form : 

n t = {{x, y), -h + b(x) <y< rj(t, x)} . 

For the sake of simplicity, we assume here that b S H /fc '°°(R d ), k being as large as 
needed, where we recall that W k '°°(R d ) = {u / d™u G L°°(R d ),0 < m < k}. In order 
to avoid some special physical cases such as the presence of islands or beaches, we set 
a condition of minimal water depth : there exists a strictly positive constant h min such 
that 

ri(t, x) + h - b(x) > h min , (t,x)eRxR. (0.1) 
We introduce the following dimcnsionless parameters : 

a hi (3 

ho \ z h 

where a denotes the typical amplitude of the waves, ho the mean depth, A the typical 
wavelength and [3o the typical amplitude of the bottom topography. In the present long 
waves regime, one has e <C 1, /x <SC 1 and e ~ /i, i.e. the Stokes number S = e/fj, is of 
order O(l). Moreover, one has (3 = Bo£ with Bo = 0(1) since we focus here on slightly 
varying bottoms, i.e. bathymetries of small amplitude. We do not make any mild-slope 
assumption on the spatial variations of the bottom, i.e. the bathymetry can vary rapidly. 
For the sake of simplicity, we take the Stokes number S and the constant Bo equal to 
one, which implies that we have fj, = (3 = e. This choice only lightens the writings by 
suppressing the constants S and -Bo from the equations, and has no influence on the 
following results. 

In [16j . the author justified a whole class of symmetric Boussinesq systems as being 
asymptotic models to the water waves problem for slightly varying bottoms in 2-D and 
3-D, assuming the existence of the water waves solutions on a large time scale. This 
assumption has been recently proved by Alvarez-Lannes in [2], where the authors sys- 
tematically justified the main asymptotics models used in coastal oceanography. The 
justification of these symmetric Boussinesq models is hence complete. All the details on 
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the construction and justification of these models can be found in [16) : the expression of 
such symmetric systems in 1-D surface is as follows : 



(£) 



rl 3 1 

(1 - ea 2 dl)d t v + d x r/ + e -r)d x r] + -vd x v - -bd x ij + a^r] 



(1 - ea^)d t ri + d x v + e -d x ((// - b)v) + a^v 



0. 



with 



a i 



-Ai- 



1 



a 2 = (Ai-1)- 



1 



2 ' v 1 ' 2 ' 

x < 02 ^ h \ u 02 1 ^ 

a 3 = A 2 (y - -) ; a 4 = (1 - A 2 )(y - g) , 

and (0, Ai, A2) being chosen in [0, 1] xR 2 such that a% — 03, a 2 > 0, 04 > 0. The parameter 
9 determines the height y = — 1 + (9 (577 + 1 — (3b) where the above velocity v is taken. The 
parameters Ai and A2 are arbitrary parameters and have no physical meaning. 
Several different set of values for (0, A, /x) can be found such that the previous condition 
on 01, d2, as, ti4 is verified, f.e. 9 = y/2/3, Xi = X2 = 1/2 which gives a\ = a 2 = 03 = 
04 = 1/12. Throughout this paper, we denote by (S) any of this symmetric system. This 
system is the starting point of our work. 



Organization of the paper 

The paper is organized as follows. In Section I, we recover the KdV approximation 
proposed by Schneider- Wayne [SD] and Iguchi [53] as follows : we first diagonalize the 
system (S) and then look for an approximation of the solution of the obtained coupled 
system. This approximation is searched under the form of a couple of waves (Uq,Nq) 
moving in opposite directions, plus correcting terms (U±,Ni) that satisfy a sublinear 
growth condition. We show that Uq and N must satisfy a system of two uncoupled KdV 
equations, as shown by Iguchi in |26j . It is proved that this approximation is correct for 
sufficiently decaying initial data and topography. 

In Section II, we discuss the validity of this approximation - which is equivalent to check 
the sublinear growth condition on the correctors {U\, N\) - for non trivial bathymetries. 
An analysis of these terms allows us to weaken the decay assumptions on the initial data 
and the bathymetry. However, it appears clearly that some more general bathymetries 
can unvalidate the model, and two examples of bottoms are provided for which the 
approximation diverges : a bottom corresponding to a simple step and a slowly varying 
sinusoidal bottom. A topographically modified KdV approximation is then proposed by 
adding bottom-dependent correcting terms to the classical approximation. 
In Section III, all these models - the Boussinesq one, the usual uncoupled KdV one 
and the topographically modified version - are numerically integrated and compared on 
the two bathymetries introduced in Section II. The numerical schemes are presented, and 
the results show that both the Boussinesq model and the alternative KdV approximation 
successfully reproduce the expected physical phenomenons. 
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1. The classical KdV approximation 



In this section, we recover the usual uncoupled KdV approximation justified by Schnei- 
der and Wayne in [50] for flat botoms, and by Iguchi in [26] for slightly varying bottoms. 
To this end, we start with any of the symmetric Boussinesq system (E) derived in [16] 
and look for approximate solutions of the diagonalized version of (E) under the form of 
two waves moving in opposite directions. Each wave is shown to be slightly modulated 
in time with a dynamic governed by a KdV equation, while the correcting terms must 
solve an inhomogeneous transport equation. The sum of these two waves is proved to 
give an approximation of the solutions of (E) for sufficiently decaying initial data and 
topography. 

1.1. Derivation of the approximation 
Let us recall the expression of (E) : 



(£) 



rl 3 1 

(1 - ea 2 d x )d t v + d x r\ + e -nd x n + -vd x v - -bd x n + aid^n 



(1 - ea±dl)d t n + d x v + e ^d x ((rj - b)v) + a 3 dlv 



0. 



with a2 > and > 0. We point out that such a system is well posed for sufficiently 
smooth initial data and provides an approximation of the water waves problem of order 
0(e) for times of order 0(1/ e) (see [TB] and [5] for further details). 
In order to recover the KdV approximation, we first diagonalize (E) by introducing the 
following unknowns : 

U = v + rj ; N = v — t] . 

Plugging the relations v = (U + N)/2 and rj = (U — N)/2 into (E) yields the following 
coupled system (r) in terms of U and TV : 



(r) 



d t U + d x U + e 



d t N ~d x N + e 



ld x (3U 2 + N 2 + 2UN) - \bd x U - -d x b(U + N) 
8 2 4 



■^d 2 x d t (U + N)- ^d 2 x d f (U — N) + a.dlU 



= 0, 



\d x (U 2 + 3N 2 + 2UN) + \bd x N + -d x b(U + N) 
o z 4 



-^d 2 x d t (u + n) + ^d 2 x d t (u -N)- ai d 3 x N 



= 0. 



At this step, we choose to look for an approximate solution (U app , N app ) of (T) of the 
form : 



U apP (t, x) = U Q (T, x-t)+ eUi(T, t, x) 



k N apP (t, x) = N (T, x + t)+ eNi (T, t, x) , 



(1.1) 



where T is the slow time variable T = et. The choice of variables x — t for Uq and x + t 
for iVg emphasizes that we look for solutions which - at first order - propagate at speed 
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one, to the right for U and to the left for N : this comes from the first order equations 
on U and N which are respectively (dt + d x )U = and (dt — d x )N = 0. 
The use of two different time scales is useful in capturing both the short time evolution 
of the wave and the nonlinear and dispersive dynamics which occur for larger time scales. 
We complement this ansatz with two initial conditions on (U app , N app ) and a classical 
sublinear growth condition on the correctors (U±, N±) : 

U apP \ t=0 = U\ t=0 , ^Vopp|t=o ~ N\ t=0 , 



and for all T > 

Aol 1 ^" L-< M.,„j.l 



'i?L7 l^i(->*>-)li«>(ro,2blii 2 (K)l -° 



| (1.2) 
lim -|iVi(.,i,.)U-([o,r ];i 3 (K)) = ■ 

t — >+oo J 

Such a condition on the correctors is quite usual in multiscales expansions and has been 
first introduced in the context of nonlinear geometric optics by Joly, Metivier and Rauch 
in [30 . It forces the correctors eU\ and eN\ to be small on the large time scale asso- 
ciated to the KdV dynamics ; namely, it ensures that eU\ — o(l) and eN\ = o(l) in 

L°°([0, y];i 2 (R)). 

Plugging this ansatz into (r) and neglecting the terms of order 0(e) yields the following 
system : 

(d t + d x )Ui = f(T, x-t)- \d x N% - \d x (U N ) + hd x U + -U Q d x b 

o 4 z 4 

+ i NAb+ ^^, (i3) 

(d t - d x )Nx = g(T, x + t)- -d x U% - -d x (U N ) - -bd x N - -U Q d x b 
-\N d x b- a -^±dlU , 



where < 



f(T, x-t) = -d T U - ^d x U* - (oi + ^^)d 3 x U , 
g(T, x + t) = -d T N - ^d x N* + ( 0l + °±±?±) d 3 x N 



(1.4) 



N.B. : For the sake of simplicity, we have kept here the notations d x on Uo and No while 
we should have written rigorously dx_ Uo and dx + Ng with X- = x — t and X + = x + t. 



At this step, an explicit resolution of (| 1 .3[) in terms of Uq and No shows that U\ and N\ 
have the simplified form : 

U 1 (T,t,x)=tx f(T,x-t) + hl(Uo(T,x-t),No(T,x + t)) , 

Ni(T,t,x) =t x g(T,x + t) + hl(Uo(T,x-t),N (T,x + t)) , 

the complete expression of (Ui,Nx) being given at the end of this section in (| 1 . 9[) . 
Such expressions of U% and Ni include terms that grow linearly in time, which is in- 
consistent with the sublinear growth conditions p.2[) . It follows that both f(T,x — t) 
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and g(T, x — t) must be null quantities for all (T, t, x) G [0, To] 



Writing explicitly 



this result, remarking that a\ + (a 2 + a^)/2 = (ai + 03 + a 2 + a^/2 = 1/6 and plug- 
ging it into (|1.3[) . one obtains the following uncoupled KdV equations on (Uq,N ) and 
inhomogeneous transport equations on the correctors (Ui,Ni) : 



d*N = 



(1.5) 



and 



1 



1 



: d x (U N ) + -bd x U + -U Q d x b 



N d x b - 



a 2 - a 4 



%N , 



1 



(1.6) 



(dt-d x )N! = --d x UZ 
8 1 



-9 x (C/ iVo) 



-N d x b - 



a 2 — a 4 



dtUn 



Finally, we can construct an approximation of the solutions of the initial system (S) in a 
natural way : let '7s)o<e<e be a family of solutions of (E) with initial data (t>o,??o)- 
One defines uq = v + 770 and no = t>o — 770, and let (Uo,Nq) be the solutions of the 
uncoupled KdV equations (Exdv) with initial data (ito,?io)- I n the end, the uncoupled 
KdV approximation of the solutions of (E) is given by : 



(M) 



v KdV 



U a +Nn 



Uo-N, 



VKdV — 







(1.7) 



This approximation corresponds to the one proposed by Schneider- Wayne [50] for flat 
bottoms and by Iguchi [26J for uneven bottoms. 

Remark 1.1 Iguchi derived his KdV approximation in the framework of capillary- gravity 
waves, and not gravity waves as in this paper. However, we can see in [261 that the only 
impact of this capillarity- gravity waves approach lies in the - constant - coefficient in 
front of the dispersive terms. The comparison between this work and the present one is 
hence far from being inappropriate. 



1.2. Validity of the approximation for sufficiently decaying topographies 
Recalling that (U app , N app ) is given by (|1.1|) . we introduce the following quantities 



v„ 



U a pp ~f~ N a pp E U a pp N a pp 



Vapp - o ' (.-^J 



app 2 ' 1 app 

in order to prove the following proposition. 

Proposition 1.1 Let s > 2, a > s + 5, (v 0) r) ) £ H a (R) and b £ W 1,00 (M). There 

T 

exists To > and a unique family T]£)o<e<e S L°°([0, — ];iJ CT (IR)) of solutions of 

(E) with initial data (vq, tjq). We define (uo, uq) = (v$ + r/o , vq — r/o) . Then there exists a 
unique solution (Uq,Nq) to the system (YiKdv) with initial data (uo,rio) and this solution 



is bounded in L°°([0, T ]; H° 

T 

Moreover, we have the following error estimate for all t € [0, — ] 

£ 



L°°([0,t];if s (R)) 



< C(l + |(t r i,A r i)| L oo ([0iTo]x[0!t] . ff s+ 3(R)) )e 2 i. 



where (uJL,?7<JL) «s defined in U.8]) . 



PROOF. The result on the system (£kw) is a very classical result on the KdV equation 
that has been established f.e. by Bona and Smith ([13], 1975) and we omit the proof here. 
The leading terms (Uq, No) and the correcting terms (f/i, N\) have been chosen such that 
(U ap p, N app ) is solution of the system (T) with a residual of order 0(e 2 ). This residual 
denoted by (e 2 i?i, e 2 i?2) can be computed explicitly and we get : 

Ri = -\d x {U U x ) - ^d x (N Ni) - jd w (U Ni) - ^d x {UxN ) + h>d x U x 
+\{U X + N x )d x b - ai d x U x + a -^±dld t U l + ^±dld t N x 



— d x d T U + — - — d x d T N 



\d x Uf + ^d x N 2 + ^(t/iiVi) 



with a similar expression holding for the residual i?2 of the second equation of (T). We 
first use {T*Kdv) and (S corr ) to express OtUq, OtNq, dtU\, dtN x in terms of spatial 
derivatives of Uq, Nq, U%, N x and then use the following standard estimates : one first 
applies the differentiel operator d x to each relation giving R\ and i?2, then multiplies 
them by respectively d x R\ and d x R,2, and finally integrates on M. Using the fact that 

T 

s > 2, this yields easily for all t E [0, — ] : 

|(-Rl,-R2)| i oc([ 0)t ] ; H=(M)) - ^V 1 + K^l'-^l)lL«>([0,ro]x[0,t];H 8 + 3 (M))) ' 

with C depending only on \(U , N )\L^([o,T a };H 3+5 (R)) and \b\ w i,^^ R y Inverting the diag- 
onalization by plugging the relations U app = v app + i] app and N app = v app - i] app into (T), 
we easily deduce that (v app , r] app ) is solution of the system (£) with a residual bounded 
by Ce 2 (l + ^Vi)lioc([ t ]x[o t]-i/ s + 3 (R)))- Standard energy estimates - as described 
above - applied on the symmetric Boussinesq system (S) yield : 

KviiVl) ~ («a^>^a|p)k~([0,t];H=(m)) < C(l + | (Ui , N X ) | Loo ([0To] x [Qt] , H S+3(R )) ) S 2 t , 

which ends the proof. 

An easy extension of this proposition is the following corollary which gives an error 
bound for the KdV approximation. 



Corollary 1.2 Under the same hypothesis as in Proposition ] we have the following 

T 

error estimate for all t £ [0, — ] : 

{vl>Vz)-(v KdV ,VKdv) L „ m] . Hsm) - C '( 1 + l( C7 l' iV l)li-([0,T ]x[0,t];J?»+3(K))) e ( 1 + £ *)> 
where {vKdV'WKdv) * s the uncoupled KdV approximation defined in JJ.Tp . 



PROOF. One has : 

( W S; Vt.) ~ ( v KdVi VKdv) = ( W £i ^s) — ( v appi Vapp) + ( v appi Vapp) 

~( v KdV> VKdv) 

Using this relation and the error estimate coming from Proposition 11.11 yields the result. 



This corollary clearly states that the validity of the uncoupled KdV approximation only 

depends on the control of the correcting terms {U\,N\) in L°° ([0, T ] x [0, t]; H s+3 (R)^j 

norm on the large time scale 1/e. From now on, these correctors become the center of 
our analysis. 

As we saw earlier, the inhomogeneous transport equations that govern the evolution of the 
correctors (U±,Ni) can be solved explicitly in terms of Uq and Nq. Using the fact that the 

solution of the equation (d t + d x )u = / is given by u(t, x) = f(x — t + s)ds + u(0, x — t) 

Jo 

and that Ux(t = 0) = 0, we thus get the following expression : 



IMT, t, x) = -1(7V 2 (T, x + t)- N 2 (T, x-t)) + ^±(d 2 x N (T, x + t) 
lb 4 

-d 2 N (T, x -t))-^U Q (T,x- t)(N (T, x + t)-N {T,x-t)) 
+hj {T,x- t){b{x) -b(x-t))- ^d x U (T,x- t) 

1 r* 

N (T,x - t + 2s)ds + -d x U (T,x - t) / b(x - t + s)ds 

2 In 



1 



+ - / d x b(x-t + s)N (T,x-t + 2s)ds , (1.9) 

and a similar expression holds for N\ : 

Nt(T, t, x) = -1([/ 2 (T, x - t) - U 2 (T, x + t)) - ^-^(d 2 x U (T, x - t) 
16 4 

~d 2 x U (T, x + t)) - ^N (T, x + t)(U (T, x-t)-U (T,x + t)) 

~N (T, x + t)(b{x) - b{x + t)) - ^d x N {T, x + t) 

1 f* 

U (T,x + t~2s)ds d x N (T,x + t) / b(x + t - s)ds 

o 2 J 

I 



d x b(x + t- s)U (T,x + t- 2s)ds , (1.10) 



We here only deal with the case of Ui since all the method can easily be adapted to the 
case of N%. 
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The corrector U± is analysed in the following way : let Tq > 0, s > 2, <j > s + 5 and 
(uo, no) in H a (M.). We know that the solutions (Uq, No) of the problem (YjKdv) with initial 
data (uq, rig) are bounded in L°°([0, T ]; iJ CT (R)) 2 . We suppose here that the bottom 
topography b is bounded in Under these circumstances, it clearly appears that 

the first four terms of the expression of U\ are bounded in L°°([0, T ] x [0, t]; H S+3 (R)). 
Only the last four terms can be problematic and deserve a precise treatment. 
If (Uq,Nq) and b come with a sufficient decay rate at infinity, we can straightforwardly 
control these terms. To this end and following [50 , we introduce the following weighted 
Sobolev space H s ' a for all s G N and a e K : 



H>>° = {u G H"(R) I \u\%., a + x 2 )"/ 2 



< oo 

L 2 (R) 



fe=0 

We can now state our first theorem on the validity of the approximation for sufficiently 
decaying initial data and bottom topography. 

Theorem 1.3 Lets >2,a> s + 5, {v ,Vo) € H a ^(R) 2 and be H S+4,1 (M). There exists 

To 

To > and a unique family (w^, ?7s)o<e<e ^ L°°([0, — ];_ff cr (R)) of solutions o/(E) with 
initial data (vo,r)o). We define (uo,no) = (vq + r/o, i»o — Vo)- Then the solution (J7 0) JVo) of 
the system (TiKdv) with initial data (uo,no) is bounded in L°°([0,To];H a ' (M)). More- 

T 

over, we have the following error estimate for all t £ [0, — ] : 

Vs) - (vKdV'VKdv) T ,, n . „ „,„ < Ce(l + st) , 

L°°([0,t];J? s (K)) 

where (v Kdv ,r] Kdv ) is the uncoupled KdV approximation defined in \1. 9j ). 



PROOF. We know from [33] and [50] that the KdV equation propagates the regularity 
of initial data taken in weighted Sobolev spaces and we omit the proof here. The end 
of the proof is devoted to the estimate of \(U\, -/Vi)| i oo([ ,T ]x[o,t];ff s + 3 (R]) 2 - The work of 
Lanncs in [39] is here very useful to control this quantity. Indeed, using the equations 
(T, corr ), the fact that Uq(T, .), No(T, .) and b are bounded in H a,1 (R), and Proposition 
3.5 of [39], one finally obtains the estimate : 

\(Ui, Ni)\ L a ^o,T ]x[0,t];H 3 + 3 m) 2 - CdH-f/^fTt)' (^0' No) L» ([0,t];H< r > 1 (R)) 2 ) ■ 

Pluging this last estimate into the result of Proposition 1.1 ends the proof. 



Remark 1.4 As specified in U^f , this approximation diverges on a large time scale in 
the periodic framework unless we specify a zero mass assumption on the initial data uq 
and no- This drawback is dealt with at the end of the next section. Until then, the results 
provided can be extended to the periodic framework with this zero mass hypothesis. 
Remark 1.5 It is worth pointing out that the validity of the uncoupled KdV approxi- 
mation for the Boussinesq system (£) is enough to demonstrate its validity regarding the 
water waves problem. Indeed, we can deduce from J16f that the error estimate between 
the solutions (wsj 7 ?s) °f (^) an( ^ ^ e solutions of the water waves problem is of order 
0(e(l + st)). An error estimate between the solutions of the water waves problem and 
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the KdV approximation (vKdv^Kdv) can thus ^ e immediately deduced from the results 
of this paper. 



2. A topographically modified KdV approximation 



In this section we discuss the validity of the previously derived uncoupled KdV ap- 
proximation on a large time scale for different bottom topographies. We demonstrate its 
validity for less restrictive bottoms, but provide two examples of simple bottoms for which 
the approximation diverges. A new approximation that takes the bottom into account is 
finally derived. 



2.1. Discussion on the validity of the approximation 



Starting from the previous theorem, it is worth wondering if this one holds for less 
restrictive initial data and bottoms, i.e. without any condition of a sufficient decay rate 
at infinity. In this view, we focus in a more general way on the last three terms of U\ 
by supposing that (uo,n ) is bounded in L°°([0, t]; H <T (R)) 2 , which is propagated by the 
KdV equation on (Uq, N ) (see [M]). Using Cauchy-Schwarz inequality on the first two 
terms and Proposition 3.2 of [39] on the last term, we can write the following controls 
T 

for all t £ [0, — ], s > 2 and a > s + 5: 



d x U (T,. 
d x U (T,. 



N Q (T,.-t + 2s)ds 



b(.-t + s)ds 



< CtVt, 

< °2 Mlhm.) . 



f d x b(.-t + s)N Q (T,.-t + 2s)ds < C 3 \d x b\ H ., R) Vt , 

JO if s (K) 

where the constants Ci, 62,63 depend exclusively on | (t^o, No)\l°° ([o,t}-,H° W) 2 ■ 
These preliminary estimates are at the heart of the proof of the following theorem. 

Theorem 2.1 Let s > 2, a > s + 5, (u ,?7o) e H a (R) 2 , b e H S+A (R). There exists 

T 

To > and a unique family (w^, ?7s)o<e<eo ^ L°°([0, — ];_ff cr (R)) of solutions o/(E) with 
initial data («o, rjo) . We define (uq, uq) = (vq + 770, v — 770) • Then the solution (Uq, Nq) of 
the system (EkcIv) with initial data (uo, uq) is bounded in L°°([0, To]; H a (M)). Moreover, 

T 

we have the following error estimate for all t £ [0, — ] : 



where (v KdVl r] Kdv ) are as defined in \1. 7| ). 



< CeVi(l + et) 



PROOF. Using the three previous inequalities, one obtains : 

\(Ui, -/Vl)|L°°([0,T o ]x[0,t];i? s + 3 (K)) < CVt 
11 



where C = C(|&|#s+4( ffi ), (Uq, -/Vo)L~([o,t];ff<'(R)) 2 )- The final result follows from Corollary 
PI 



This theorem proves that the approximation is less precise on a large time scale if we 
remove the assumption of a sufficient decay rate at infinity. And yet, it is worth pointing 
out that the regularity imposed on b in this theorem excludes many physical cases of 
interest. We focus from now on two simple examples of bottoms which do not fall into 
the scope of Theorem 12. II : a regular step, and a slowly varying sinusoidal bottom. Our 
goal is to emphasize the fact that the approximation (v Kdv ,r) Kdv ) diverges from the 
exact solution (u^,ry£) in these two simple cases. To deal with such bathymetries, a 
topographically modified KdV approximation is derived at the end of the section. 

In order to simplify the analysis, we only consider the approximation corresponding to 
ai = 1/6, a,2 = 0, a4 = which is obtained for 9 = \/2/3, Ai = 1, A2 = 1, and the case of 
a wave propagating to the right. This last condition is realized by taking no = 0, which 
implies that N = N = 0. 



2.1.1. The case of a step 

We consider here a bottom whose shape corresponds to a regular step. The interest of 
such an example is that in this case, b L 2 (M). 
The bottom is defined as follows : 



, Vx < 



.4 



m = { 2 ( 1 



(7(^-5))) > Vx€[0,l], 



(2.1) 



A , Vx > I 



For a right going wave, the system (Sxdv) is reduced to the simple KdV equation : 

d T U + ^d x Ui + ^dlU = , 

and we chose the initial condition uq such that the solution of this equation is a positive 
soliton which propagates to the right. 

We write the explicit expression of the corrector U\ when Nq = : 



Ui{t,x) = hj {x-t){b{x)-b{x-t)) + -d x U Q {x-t) 



b(x — t + s)ds 



In this expression, the only possibly secularly growing term is d x lIo(T,x — t) / b{x — 

Jo 

t + s)ds. The time evolution in amplitude of this term is obviously led by the evolu- 
tion of / b(x — t + s) for all x G M. When the bottom is a step as defined in (|2.1[) . 
Jo 

this integral essentially grows linearly in time. We now prove that because of this, 
|t r i|i°°([o,T ]x[o.t];ff s + 3 (R)) grows linearly in time. Let s > 2 and er > s + 5. Starting 

from the expression of U\, we get for all t <G [0, — ] the following estimates : 

e 
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|*7i(T,V)|h«+3(r)> -d x U (T, ■ - t) / b(--t + s)ds 



ff 3 + 3( 



c 



with C = 



~f/ (T, ■ - *)(&(■) - &(• - t))|^ +j < i|ft|L-|^0|z-([0,t];ff.+3) = Co, 



|f/ 1 (T ) i,-)| H »+3 ( E) > ^|fl*Ub(T, • -*) / ''(■- /-:-«)</>■! l^,.., - r„. 



\d x U (T,x-t)\ 2 / 6(a;-t + s)ds 



— Cn 



since / 6(a; — t + s)ds — , Vx < 
Jo 



> 



= -At , 



\d x U (T,x-t)\ 2 / 6(s)ds 



l+i 



dx — Cq , 



\d x U (T,x-t)\ 2 dx-C , 



r 

since / b(s)ds = At , Vx > 2 + £ 
J x -t 



1 



which implies that 



l^llL°°([0,T o ]x[0,t];ff s + 3 (R)) - Cl* Co i 



where the last positive constant C\ only depends on |(? x t/o|i/ 2 ( 



(2.2) 



This linear growth of |£^i|z,°°([o,T ]x[o,t];-ff s + 3 (R)) i s sharp since it follows from the explicit 
expression of U± that this growth is at most linear. Furthermore, we recall that 



{viiVi) ~ ( v Kdv, View) = ~ (w a p P , V<£p) + £ (— 



m + m Ui- Ni 



(2.3) 



Using this relation, (|2.2p and Proposition ll.il we get that there exists two constants C2 

T 

and C3 and a time Ti independent of e such that Vt G pi, — ], 

e 



{vi,r)D- (VKdViVKdv) 



L°°([0,t];ff s 



> |C2(l + t)e-C 3 (l + t)£ 2 d 



We finally deduce that there exists two constants C and C" such that Vi S [Ti, — 1, 

e 

This proves that in this case, the error is of order 0(1) on times of order 0(1/ e), and 
the usual KdV approximation is not valid for such a topography. 
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2.1.2. The case of a sinusoidal bottom 

We consider here a bottom denned as follows 



b(x) = Asm(sx) , VieR. (2.4) 
where Uq is again a soliton propagating to the right. 

We mention that such a type of periodic bottom varying on a slow spatial scale has been 
studied in [TS] by Craig-Guyenne-Nicholls-Sulem, with the difference that the authors 
authorized the bottom to vary also on a small spatial scale. 

Again, the amplitude of the term d x Uo(T, x — t) I b(x—t+s)ds evolves in time according 
t 

to / b(x — t + s)ds. Let us have a look at this quantity for all x € K and t > : 
Jo 



b(x — t + s)ds = / sin(sx)ds 

Jx-t 

= [cos(ex) — cos(e(x — t))] 

e 

2A , ( t \ . (et 
= — sin e{x - -) sin — 



We can see that the amplitude of this term is of order 0(l/e). We now demonstrate that 
it is also the case for the corrector U\ : 



\Ui{T,t 



' 'J\H'+ 3 ( 



> 



\d x U {T, . - 1) / b(- - t + s)ds 
* Jo 



Co : 



H"+ 3 < 



> 



d x U (T,--t) I b(--t + s)ds 
o 



Co ) 



i v r \ d * uo{ ?> x - ^ 2 sin2 ( £{x - sin2 (?) dx - c ° ' 



Ism 2 . 2 



| ( 7,T'„(r..r -t)| 2 sin 2 (e(x-|))«/.r. r„ 



(2.5) 



At this point, we remark that 

0<J \d x U Q (T,x-t)\ 2 sin 2 (e(x-^dx< j \d x U {T, x ~ t)\ 2 dx , 
We hence deduce that for all t > there exists a(t) £ M such that 

\d x U (T,x-t)\ 2 sm 2 ( e(x - t)) dx = sin 2 (a(t)) f \d x U a {T, x - t)\ 2 dx . 



Pluging this one into (|2.5|) leads to 

A 



\U x {T,t,-)\ He+H 



> 



snr f y ) sin 2 (cm/ i 1 



d x U (T,-) 



Co 
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which finally implies that there exists a constant C\ such that 



|^iU~([o,T ]x[o,t];H*+ 3 (m)) > - Co ■ 



Using this result and the same technique as in the previous example leads to the same 
conclusion : the uncoupled KdV approximation diverges on a large time scale in this case 
too. 



2.2. A topographically modified approximation 

Both examples clearly show the invalidity of the approximation on a large time scale if 
we consider general bottoms topographies b which do not have specific decay properties 
at infinity. Therefore, we need to modify the usual KdV approximation to be able to 
handle general bathymetries. 

All the previous analysis has shown that two terms of the r.h.s. of the explicit expression 
(|1.9[) of U% may exhibit a secular growth, and cause the approximation to diverge on a long 



1 r 1 r 

time scale: these are —d x U$(T, x — t) I b(x — t + s)ds and - / d x b(x — t + s)N (T,x — t + 

2 Jo ^ J 

2s)ds. As far as the expression (|1.10| of Ni is concerned, the same possibly problematic 

terms are -^-d x No(T, x + 1) 

I b(x+t— s)ds and — — / d x b(x + t — s)Uo(T, x + t — 2s)ds. The idea is as follows: rather 
Jo 4 J 

than treating these terms as correcting terms - which invalidates the approximation for 
general bathymetries - we can include them with the leading order one terms Uq and N 
in the final approximation. 

This idea leads us to propose the following topographically modified KdV approximation 
which is an alternative version of (M) : 



15 



e,b 
v KdV 



U Q + N a 



(M b ) { 



e,b 
^KdV 



Uo-N 



d x U {T,x-t) / b(x-t + s)ds 
t 

-d x N (T,x + t) / b(x + t-s)ds 

t 
+- / d x b(x-t + s)N (T,x-t + 2s)ds 
2 Jo 

i r* 

— — / d x b(x + t - s)U (T, x + t- 2s)ds 
+-U (T,x-t) (b{x)-b(x-t)) 
+ ^N (T,x + t) (b{x + t) - b(xj) 



d x U {T,x-t) / b(x-t + s)ds 



-d x N (T,x + t) / b(x + t-s)ds 

+ o / d x b(x-t + s)N (T,x-t + 2s)ds 



1 



1. 



d x b(x + t - s)U (T, x + t- 2s)ds 



+-U (T,x-t) {b(x)-b(x-t)) 
+ ^N {T,x + t) (b{x + t) - b(xj) 



(2.6) 



where Uq and Nq are still solutions of the system (H K e dv ), and where all the possibly 
secularly growing terms of (|1.9j) and (|1.10|) have been included with the order one terms. 
The physical role of each of these additionnal topography-dependent terms is discussed 
in the last section, where we validate this model numerically on the previous examples 
of a step and a sinusoidal bathymetry. 

Remark 2.2 We have here also included the terms Uo(T,x — t)(b(x) — b(x — t)) and 
Nq(T, x + t) (b(x + 1) — b{x)) even if these terms remain bounded indepently of e for all 
time. The reason of this choice is that we are interested in their physical meaning. Indeed, 
we further see - in the last section - that they are responsible for the reproduction of the 
phenomenon of shoaling. We hence decided to include these terms in the approximation. 

The main advantage of this modification relies in the following remark : now that the 
bottom terms have been included with the leading order terms in the approximation, we 
can easily see that the correcting terms U\ and iVi solve a different equation. Indeed, the 
equations on U\ and N\ become : 

{dt + d x )Ui = -ld x N* - ]d x (U N ) + ^-^i d*N , 

(Kir) { 11 a -a (2J) 

ft - d x )N x = ~-d x U 2 - -d x (U N ) - d x U Q . 
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It is clear here that all the possibly secularly growing terms of the correctors (Ui,Nx) 
have been removed. 

Remark 2.3 Numerically speaking, this modified version is quite interesting since the 
topographical terms are computed explicitly from the solution of the KdV equations. We 
thus expect the numerical simulation of this model to be faster than the one of the sym- 
metric Boussinesq model (£). This point is checked in the last section. 

In the periodic framework, we saw that the usual approximation is not valid on a large 

time scale because of the linear growth in time of the term d x Uo(T, x — t) / Nq(T, x — 

Jo 

t + 2s)ds in U±, unless we specify a zero mass assumption on the initial data uq and 
no- Once more, we can propose a valid approximation just by including this term in the 
order one terms of the ansatz. We conclude this section with the proposition of a new 
approximation in the periodic framework : 



v 



e,bper 

KdV ~ 



(M pe b r ) 



e.b per 
" KdV 



e,b _ 

v KdV o 



d x U (T,x-t) I N (T,x-t + 2s)ds 
o 



+d x N (T,x + t) / U (T,x + t-2s)ds 



e.b 

VKdV g 



(2.8) 



d x U (T,x-t) I N (T,x-t + 2s)ds 
o 



-d x N {T,x + t) / U (T,x + t-2s)ds 
Jo 



To end this section, we would like to add a few words on the works - referenced in 
the introduction - 48 32 28 43 45] and [55]. These papers are mainly devoted to the 
derivation of KdV equations in a variable medium, but differ from the present work in 
the way the topography is accounted for. Indeed, the bathymetry does not show up as 
additionnal terms computed from the solutions of the uncoupled KdV equations like here, 
but as variable topography-dependent coefficients in the KdV equations themselves. A 
good example of this result is the coupled KdV model obtained by Iguchi in [26] . Their 
approach is thus very different from ours, and lead to additionnal difficulties such as 
the well-posedness of these topography-dependent KdV-like equations. However, Iguchi 
[21)] has been able to show that his coupled KdV approximation is well-posed in some 
Sobolev spaces and valid - in the meaning of the precision of the approximation given by 
the solutions - on a long time scale, as long as - among other conditions - the topography b 
belongs to W k '°°(R) for k large enough. This model is hence quite pertinent in the present 
context of realistic bathymetries and offers a good alternative to our topographically 
modified KdV approximation. 

Nevertheless, it is worth pointing out that there actually exists a way to recover variable- 
coefficient KdV equations as in [48 32 28 43 45 126] . simply by modifying the ansatz (jl.ip 
as follows : 
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U app (t, x) = (1 + ~ b) U (T, x - (1 - | 6)t) + t, a;) , 

(2.9) 

N app (t, x) = (1 + | 6) iV (t, x + (1 - | 6)t) + eJVi(T, t, a:) , 

Plugging this ansatz into (T), neglecting all the 0(e 2 ) terms but the e 2 bd x U\ and e 2 bd x Ni 
ones, and then proceeding as in section 1 leads finally to the same KdV equations on Uq 
and No, but to a new system on the correcting terms : 

- | b)d x )Nt - -i&C/ 2 - ^(tWo) - \uod x b - figff . 

(2.10) 




Finally, defining 



i(t,x) = (i + ^b)Uo(r,x-(i-^b)t) , 

i(t,x) = {l+~b)No(T,x + (l-~b)t) , 



yields the final set of uncoupled KdV equations 
d t u + d x u + e 

dtn — d x n + e 
and the final KdV approximation 



—ud x u H — d x u bd x u d x b u 

(Sj&v) \ ^ \\ \ \ \ \ (2.H) 

—nd x n——din+—bd x n+—d x bn 
4 6 2 4 . 








, j = ^- + g I d x b(x-Pat-s))N (T,x + s-P b s (t-s))ds , 

^ = ^ " I / ^ + n e (i - »))^o(T, x - s + i?(t - a))da , 



(2.12) 

with P b e = 1 — — b, and where we need to include the possibly secularly growing terms 

in the approximation. This time, the derived KdV equations are topography-dependent 
and look very similar to the one derived in f.e. [48 32 28 . This alternative formulation 
remains uncoupled and allow us to establish a link between our approach and the older 
works on the KdV approximation over uneven bottoms. 



3. Numerical comparison of the models 

This section is devoted to the numerical comparison of the different models involved 
in this article. We compare here three models : the symmetric Boussinesq system (£) 
coming from [16] . the usual uncoupled KdV approximation justified by Schneider- Wayne 
([50j. flat bottoms) and Iguchi ([26], uneven bottoms), and finally the topographically 
modified KdV approximation. The aim is here to compare these three models for two 
non trivial examples of topography : a step and a slowly varying sinusoidal bottom. 
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3.1. Numerical schemes 



Our goal is to compare three models, the symmetric Boussinesq one, the usual KdV 
approximation (M.) and its topographically modified version (Mb)- The comparison is 
made for a solitary wave propagating to the right above two topographies : a step and 
a slowly varying sinusoidal bottom. We use for the Boussinesq system (E) and the KdV 
equations (T,Kdv) a Crank-Nicholson scheme combined with a relaxation method coming 
from Besse-Bruneau in [8] and justified by Besse in [6]. This type of scheme is of order 
two in space and time, which is appropriate for our purpose. 



3.1.1. Numerical scheme for the KdV approximation 

Due to the identical structure of the two KdV equations of (£#-<iv), we only present the 
numerical scheme for the first equation. Defining u(t, x) ~ Uo(T, x—t), we can reformulate 
this equation as follows 



d t u + d x u + e 



1 



3 , , 

-ud x u + -<9> 



= 



(3.1) 



We use a Crank-Nicholson scheme and the relaxation method introduced by Besse- 
Bruneau in [8] and justified by Besse in |6J which replace the costly numerical treat- 
ment of the nonlinear term by a predictive step. This provides us with the following 
semi-discretized in time equation : 



,n+l 



(It 



,n+l 



if „_ Li „ fu n+1 



+ (1-a)' 



.71+1 



1 n n+1 ~ 
-d 3 - 

6 x 2 



2 

-u r - 



0, 



where the predictive term u n+2 is defined as follows 



,n+4 



1-1/2 



The discretization of the nonlinear term ud x u here takes advantage of the two possible 



discretizations 



n+ 



2 d x 



71+1 



and 



a 



71+1 



u 



a, 



u 



by introducing a parameter 



2 J 2 

a £ [0, 1] and taking a convex combination of these possibilities. Keeping in mind that 
we want to preserve the semi-discrete L 2 norm, an easy integration by parts gives us the 
appropriate value a = 2/3. We then choose the spatial discretization so that the discrete 
L 2 norm is preserved by the scheme, which gives the final discretization of (|3.1[) : 



St 



D x - 




(3.2) 



Dx- 



= o 



where the matrix D\ and D% are to the classical centered discretizations of the derivatives 
d x and 9 . 
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Once the solution of the KdV equation is obtained, the KdV approximation is built 
thanks to the relations (|2.6[) . by computing the integrals corresponding the correcting 
terms with the composite trapezoidal quadrature rule. Here we have made the choice to 
save the solution u of the KdV equation, and then build the approximate solution from it, 
which can be computationnaly costly in general, especially in 3-D. This is not problematic 
here since we work in 2-D. However, a optimised approach would be to see the equations 
(|2.6[) as an easy to discretize two dimensionnal pseudo-differential operator acting on 
the vector (Uo, No). Multiplying the previous scheme with the discretized version of this 
operator would lead a new system to solve at each time step, which would give us directly 
the KdV approximation. 



3.1.2. Numerical scheme for the Boussinesq system 

As far as the discretization of the Boussinesq system (S) is concerned, we consider the 
same ideas. Using a Crank-Nicholson scheme and the same relaxation method, we aim 
here at preserving the specific norm \{v,rf)\ H i = \v\^2 + MI2 + £a2\d x v\x J 2 + eai\drj\ L 2 ■ 
This quantity is indeed conserved by (E) (see [12] for more details). To this end, the 
nonlinear terms vd x v, rjd x r), r\d x v and vd x r\ are discretized in order to preserve both this 
specific discrete norm and their symmetric structure. Remarking that the equalities 



(vd x v, v) L 2 = ; (r)d x r}, v) L 2 + (r)d x v, ij) L 2 + (vd x r], -q)^ = , 



hold for (E) and using the same kind of method as for the KdV equation leads to the 
following semi-discretization of the nonlinear terms : 



< 



vd x v(nSt) s 


2 n+ i. ,v n+1 ■+ 
* r + *d x ( 2 


— )^ 


1 v n+1 -+ 
3 2 




r]d x r](n5t) p 


1^ 2d *( 2 


■ „y 


1 n n+1 4 
^3 2 




rjd x v{n5t) s 


1" ^ 2 


- vn )- 


1 TJ n+1 -i 

3 2 


vn d x v n+ i 


vd x i](nSt) s 


2 n+ i_ ,< +1 4 
V ' 4( 2 




1 v n+1 -i 
h 3 2 


' vn d xV n+ i 



We then choose the spatial discretization so that the discrete H e norm is conserved, and 
these ruminations yield this final scheme : 
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(I-ea 2 D 2 )- 



+e 



Mi 



St 



v n+l + v n 



e n n+1 + n n 

V-\b)dJ—±1- 



M-, 



+01 D 3 



= 



(7 - ea±D 2 ) 



r) n+1 - rj n 



+e 



M 



St 



v n+l + v n 



+ [(I-- 2 B) Dl 2 



Ma 



+01 £>3 



w n+l + y n 



(3.3) 



where the matrix (Mj)i<j<4 are as follows : 



Mi 



M 



77"+ 1 + 77" 



M: 



v n+l + v n 



M 



^n+l + 



2 (t, < + 



i+1 + u i-l 



) U>i 



+ 



+ 7J? 



k +3 1 o 



r? «+ 1 + 



+ 



5 (^*) 



v? +1 + v? 



iruti 2 +vti 2 f D v n+1 +v n 



I „ 1 1 



™+5 , U i+1 



{ "To . 



— ) ( D i 



v n+l + y n 



n+1 + n 



The matrix £>i and U3 are as defined in the KdV scheme, and the matrix D 2 is the 
classical centered discretization of the derivative dt. 



3.1.3. Initial data 

Let us now talk about the initialization of the two schemes. First, all the prevision terms 
are initialized with a simple explicit integration of the equations on a half-step in time. 
Then, the initial conditions are chosen such that the simulated wave is unidirectional and 
propagating to the right. To this end, we first take the initial data of the second KdV 
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equation to be zero. Then the system (T.Kdv) reduces to the equation (|3 . 1 1) for which we 
know the existence of solitary waves expressed as follows : 

^ ' ^ cosh 2 (k{x — ct + Z)) ^ 

• i 1 ea 1 3a 

with c = H , k = \ — and a, I bane arbitrary. 

4 V 8 

It is hence natural to specify the initial condition for the KdV equation (|3.ljl as follows : 

(X 

<t = ^=M*) = cosh2{k{x + d) y (3-5) 

Finally, and because of the way the KdV approximation was constructed from the Boussi- 
nesq model, we specify the initial conditions for this latter as follows : 

v(t = 0, x) = 7](t — 0, x) — — uq(x) . 

3.1.4. Validation of the numerical method 

With the initial data (|3.5|) , the KdV scheme is expected to propagate the corresponding 
solitary wave to the right, without any deformation for all time. In order to validate this 
scheme, the numerical results obtained with the initial data Q3.5P have been compared 
with the analytical solution (|3.4[) . The following relative errors on the free surface have 
been computed in the L°° norm for several values of epsilon and for computation times 
T = 1/e : 
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0.03 


1.5546. 10 -3 
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10 


80 


0.04 


0.04 


1.3717.10~ 3 
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80 


0.05 


0.05 


1.0534.10~ 3 



where L is the length of the computational domain and Sx, St are respectively the spatial 
and time discretization steps. These results allow to validate the scheme proposed for 
the KdV equations. 



3.2. Numerical results and comments 



3.2.1. Numerical results 

As specified in [16] , the choice of the parameters ai , <Z2, 04 is very interesting in a 
numerical point of view. Indeed, the parameter a\ controls the presence of the dispersive 
terms d^.v and d^.rj whereas the parameters a-i and 04 correspond to the terms d^dtv and 
d^dtT]. These last terms have the main advantage of being regularizing terms analytically 
and numerically speaking, they smooth in some way the solution because they provide a 
control of the quantities d x v and d x r\ in the L 2 norm. We decided to use here the system 
(E) corresponding to a\ = 02 = 04 = 1/12 because it is likely to provide the better 
results. 

All the forthcoming results are expressed in non-dimensionalized variables. We recall 
that both the free surface and the bottom are of size e : y — erj for the free surface and 
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y = — 1 + eb for the bottom. However, in order to get clear and readable results, we have 
plotted a rescaled free surface y = 77 and a rescale bottom y = —1 + 6. A quick word 
on the duration T of the simulations : the previous secion and [16] provide us with a 
justification of the models on large time scales of order 0(l/e). We have decided - only in 
the first example of the step - to overtake this large time scale and simulate the models 
on the very large time T = 1/e 3 ^ 2 , in order to see if the model remains stable on such 
time scales. 

The three models have been tested on two different examples of bottom. The first one 
correspond to a step at the bottom, defined similarly to [22] by 



b(x) 



, Vx G 



' 2 2 



00 {-, , • 1 L 
y (l + sm(-(x--) 



3 L 
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(3.6) 
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2 + 2' L 



where 0o is an arbitrary constant of order O(l) and L is the length of the computation 
domain. The second example corresponds to a slowly varying sinusoidal bottom, defined 
as : 



b(x) = 6 sin [ - + — x 



(3.7) 



where I is defined by I — 



1 + ea/4 



and a is the amplitude of the initial data defined in 



The following results show the snapshots of the simulations at different times - so that 
the time evolution is relatively visible - and the evolution of the relative L°° error be- 
tween the free surfaces obtained with the Boussinesq model and respectively the KdV 
approximation and the topographically modified approximation. The three models have 
been systematically plotted together in the same pictures in order to compare efficiently 
their respective behaviours. The numerical simulations have been performed for different 
values of e in the case (|3. 61) of a step : e — 0.05 and e = 0.2, which are typical values 
of the upper part of the range of validity of the long waves approximation. As far as 
the case (|3.7p is concerned, we simulated the models for the value e = 0.1. For all the 
simulations, the amplitude a of the initial free surface and the constant 0q linked to the 
bottom have been taken equal to 0.5. Here is a global tabular precising all the values of 
interest used in the simulations. 
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The figures [3. 2. 21 13.2.41 13.2.61 show the relative error between the computed free surfaces 
of the different models for each value of e and for the two cases of bottom. 
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Figure 3.2.1. Influence of the step for e = 0.05 
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Figure 3.2.2. Relative L°° error between the free surfaces for e = 0.05 
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Figure 3.2.3. Influence of the step for e = 0.2 
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Figure 3.2.4. Relative L°° error between the free surfaces for e = 0.2 
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Figure 3.2.5. Influence of a slow sinusoidal bottom for e = 0.1 
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Figure 3.2.6. Relative L°° error between the free surfaces for e = 0.1 
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Relative error between the numerical solutions ot the KdV approximation and the modified KdV approximation 



0.14 




01 23456789 10 

t 



3.2.2. Comments 

For the sake of readability, we call the computed waves as follows : B denotes the 
wave coming from the Boussinesq model (£), JCt op o denotes the wave produced by the 
topographically modified KdV approximation (Mb), an d finally JC denotes the solitary 
wave resulting from the usual KdV approximation (M). 

In the case of the step, we observe in figures 13". 2 . II and | 3 . 2 . 31 that for all tested values of 
e, both the Boussinesq model and alternative version of the KdV approximation succeed 
in reproducing the phenomenon of reflexion on the bathymetry : a smaller solitary wave 
appears on the third snapshots when the main wave goes over the step. This reflected wave 
propagates to the left at the same speed as the main wave. The classical uncoupled KdV 
model cannot - of course - reproduce this phenomenon seing that it does not depend at 
all on the bottom topography. Moreover, the Boussinesq and topographically modified 
KdV models successfully describe the following expected physical phenomenons : the 
shoaling which corresponds to the growth in amplitude of the wave after the step ; the 
deceleration of the wave after the step : the waves B and JC topo are behind the wave K, 
which propagates at a constant speed on the last snapshots of figures l3".2.1l and [3.2.3l ; and 
finally the loss of symmetry and the narrowing of the wave, which can be remarked by 
comparing - f.e. on the last snapshot of figure [372751 - the distances between several points 
of the resulting waves at different heights : the solitary wave /C propagates without any 
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deformation and remains symmetric, whereas the symmetry and width of B and JCtopo 
are modified by the step. All these phenomenons are the premisses of the process of 
wave breaking, and they are all successfully reproduced by the Boussinesq and new KdV 
models. 

A very interesting remark on our KdV model is that the role of each correcting term in 
the approximation can be intuitively identified, and these intuitions have been confirmed 
with several simulations - that are not presented here - in the case of the step. Indeed, 
it appears clearly on figures 13.2.11 and 13.2.31 that : 



- the correcting term — / d x b(x — t + s)Uo(x — t + 2s)ds is responsible for the birth and 
8 Jo 

propagation of the reflected wave, and for the very beginning of the shoaling effect, 

£ 

- the correcting term — Uq(x — t)(b(x) — b(x — t)) clearly reproduces the pursuit of the 
shoaling after the step, 



- the correcting term — d x Uo(x — t) / b(x — t + s)ds is reponsible for the deceleration 



and loss of width and symmetry of the main wave after the step. 

In the case of a slowly varying sinusoidal bottom - corresponding to figure 13.2.51 - the 
effect of the bottom is also clearly visible on B and JCtopo- To understand the speed 
variations of these waves, we have to keep in mind that the comparison is made with 
the wave JC which evolves as if the bottom was flat and located at the height y = — 1. 
Consequently, when B and JCtopo propagate above the downward part of the sinusoidal 
gap - see second snapshot of 13.2.51 - they cross two different areas : a first area - for a 
time t < T/4 - where the depth is lower than for a flat bottom located at y = —1, and 
a second area - for t < T/2 - where this is the contrary. This explains why the waves B 
and JCtopo are located for t — T/2 at the same position as JC : the waves B and JCtopo 
have speeded up over the first area in comparison with (JC), and then have decelerated 
over the second area. However, we can see - still on the second snapshot of 13.2.51 - that 
these waves are larger than (JC) at t = T/2 : this is due to the loss of amplitude of the 
waves B and JCtopo during this downward part which makes the waves be naturally wider. 
About this loss of amplitude, it is explained by the fact that the amplitude of the bottom 
decreases over the downward part of the sinusoid, which produces the inverse effect of 
the shoaling from the previous case of the step. In addition, we can see a reflected wave 
for B and JCtopo - a depression this time - which goes to the left at the same speed as the 
main wave, which corresponds to the same phenomenon of bathymetric reflexion as in the 
case of the step : the amplitude of the bottom decreases and thus produces a depression 
wave that propagates in the opposite direction. We can see that this depression wave is 
here larger than for a step because of the slow variations of the bottom. As far as the 
upward part - see last snapshot of 13.2.51 - of the sinusoid is concerned, all the previously 
described effects happen in an inverted way, and we finally recover three identical main 
waves for the three models. At this final point, the only remaining visible effects of the 
crossed topography are the reflected waves. 

As specified earlier, we decided to simulate the models on the very large time T — 

in the case of the step on figures 13.2.11 and 13.2.31 All these models have been proved to 

be valid on the time scale 0(1 /e) and it is interesting to check numerically their validity 




£ 
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- or not - on larger time scales. For a time T — 1/e 3 / 2 , we can observe on the last two 
snapshots of figures 13.2.11 and 13.2.31 that the wave ICt opo goes on growing more and more 
in amplitude, and that a depression wave deepens in front of the main wave. These effects 
are obviously not physical and can be explained by the fact that the size of the correcting 

£ f* 

term —d x Uo(x — t) / b(x — t + s)ds evolves in time like e/t as we saw in the previous 
4 Jo 

section : on a time T — 1/e 3 ^ 2 , this size become of order 0(l/y/e), which explains why 
this model diverges from the other models on this time scale. This is the main restriction 
of this model, in comparison with the Boussinesq one which seems to remain stable on 
very large time scales. An interesting perspective would be to look for higher order terms 

- like Wright in 54J - in the approximation to deal with this problem. 

To sum up, the results on these two examples of bottom show that both Boussinesq 
and topographically modified KdV models are able to reproduce the expected physical 
phenomenons : reflection, shoaling, loss of speed and symmetry. This is of course not 
the case for the usual KdV approximation which is independent from the bottom topog- 
raphy. Even if we can isolate the role of each correcting terms with our modified KdV 
aproximation, this one diverges when time goes over the theoretical limit time of validity 
T = 1/e. The Boussinesq does not have this drawback and remains stable on time scales 
of order 0(l/e 3/2 ). 
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